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BR(5 — » X s 7): First Next-to-Next-to-Leading Order Results 

Matthias Steinhauser a * 

a II. Institut fur Theoretische Physik, Universitat Hamburg, 
Lumper Chaussee 149, D-22761, Hamburg, Germany. 

In this contribution the first next-to-next-to-leading order corrections to the branching ratio BR(B — > 
X a j) are discussed. This includes the completion of the matching part and the fermionic corrections to 
the most important operator matrix elements. 



1. Motivation 



The rare decay of a B meson into a me- 
son containing a strange quark and a photon 
is very sensitive to new physics which is con- 
nected to the fact that in the Standard Model 
(SM) this process is loop-induced. Thus a pre- 
cise comparison of theoretical calculations and 
experimental measurements can lead to im- 
portant hints on the theory beyond the SM. 
Conventionally the decay is normalized to the 
semileptonic i?-meson decay. The correspond- 
ing branching ratio has meanwhile been mea- 
sured by CLEO [HI, ALEPH [E], BELLE [ 
Is3 and BABAR [ 4 which is summarized in 
Fig. H and leads to the experimental world 
average 2 [E] 

BR[B -» X sl , {Ry > m b /20)} = 
(3.42 ±0.30) x 10" 4 . 



(1) 



This result agrees very well with the SM pre- 
dictions [ElCj 

BR[B -> X s j, {E 7 > 1.6 GeV)] = 

(3.57 ±0.30) x 10" 4 , 
BR[B -» X sl , {E 1 > mb/20)] ~ 



3.70 x 10" 



(2) 



In Fig. Q] the latter is compared to the mea- 
surements. In the near future the experi- 
mental error will reduce to approximately 5% 
which reduces the current error by almost a 
factor of two. Thus also the reduction of the 
theoretical uncertainty in Eq. (0) is manda- 
tory. 



*To appear in the Proceedings of the International 
Conference and Loops and Legs in Quantum Field 
Theory 2004, Germany, 25-30 April 2004. 
2 The CLEO measurement from 1995 is not included 
in the average. 
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Figure 1. Comparison of the experimental 
measurements for BR[_B — > X s j], where a cut- 
off on the photon energy E 1 > mb/20 is ap- 
plied, with the theoretical predictions given 
by the (yellow) band. The full (horizontal) 
lines indicate the uncertainties of the com- 
bined experimental results (cf. Eq. JU). 

At next-to-leading order (NLO) a detailed 
analysis of the theoretical error has been per- 
formed in Ref . [ Ej . It has been shown that 
the largest uncertainty comes form the depen- 
dence on the charm quark mass which arises 
for the first time at this order. Including all 
next-to-next-to-leading order (NNLO) effects 
one arrives at an uncertainty of approximately 
7% which has to be compared with the 8% er- 
ror quoted in Eq. @ • Thus a complete NNLO 
calculation can significantly reduce the theo- 
retical uncertainty. 

Calculations for the process B — > X s -f are 
conveniently performed within the framework 
of an effective theory where the scales of order 
Mw and higher are integrated out. This leads 
to an effective Lagrangian of the form 



C 



off 



y (u,d,s,c,b) 



(3) 



AG 



V2 



hi 



v; s v cb c? + v*v tb ci] ZijPj 



where Gf is the Fermi constant, V stands for 
the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix, Pj are the operators, C?' are the Wilson 
coefficients and Zij the corresponding renor- 
malization constants. More details and ex- 
plicit expressions can be found in [[S] and ref- 
erences cited therein. 

In such an effective-theory framework the 
calculation splits into three steps: (i) the 
matching of the full and the effective theory 
at the high scale, (ii) the evaluation of the op- 
erator matrix elements at the low scale mj,, 
and (iii) the running from the high to the low 
scale. At NNLO step (i) requires the evalua- 
tion of two-loop diagrams in order to obtain 
the matching coefficients for the four-fermion 
operators which has been performed in Ref. [ 
[5]. The three-loop matching coefficients CV 
and Cs of the dimension-five dipole operators 
have been evaluated in Ref. [|H] which will be 
considered in more detail in Section [3 

As far as the evaluation of the NNLO op- 
erator matrix elements are concerned up to 
now only the fermionic corrections of 0(a 2 s rif) 
to Pi, P2, P-i and Pg are known [ 10 . At 
NLO these operators are numerically impor- 
tant which is a strong motivation to con- 
sider the numerical effects at NNLO. Fur- 
thermore, once the corrections proportional 
to nf are available, it is tempting to apply 
the so-called naive non-abelianization in or- 
der to estimate the complete corrections of 
order a 2 .. This is based on the observation 
that the lowest coefficient of the QCD (3 func- 
tion, (3q = 11 — 2n//3, is quite large and thus 
it is expected that the replacement of n/ by 
— 3/3o/2 may lead to a good approximation of 
the full order a 2 corrections. In Section [3] wc 
provide more details on this calculation. Fi- 
nally, Section 0] contains our conclusions. 

2. Three-loop contribution to C-j and C% 

There are several possibilities to evaluate 
the matching coefficients. For instance, it is 
possible to evaluate the Feynman diagrams 
on-shell. However, we find it more convenient 
to follow the procedure outlined in Ref. [ |H] 
where all the necessary diagrams are evalu- 
ated off-shell, after expanding them in the ex- 
ternal momenta. In intermediate steps spu- 




Figure 2. One of the 0(1O 3 ) three-loop dia- 
grams that contributes to b — ► s"f. 



rious infrared divergences appear which are 
generated by the expansion and are regulated 
dimcnsionally. They cancel out in the match- 
ing equation, i.e. in the difference between 
the full SM and the effective theory off-shell 
amplitudes. 

The scalar three-loop integrals are evalu- 
ated with the help of the package MATAD [ 
lllj designed for calculating vacuum diagrams. 
The fact that MATAD can deal with a single 
non-vanishing mass only is not an obstacle 
against taking into account the actually dif- 
ferent masses of the W boson and the top 
quark. Expansions starting from m t — My/ 
and m t 3> Myy allow us to accurately deter- 
mine the three-loop matching conditions for 
the physical values of mt and Myy. 

The matching coefficients are obtained 
from the requirement that the one-particle- 
irreducible Green functions in the full and ef- 
fective theory are equal. In the latter case 
only tree-level diagrams contribute as after ex- 
pansion in the external momenta only scale- 
less integrals appear which are set to zero in 
dimensional regularization [ I12| . However, the 
counterterm contributions that arise from the 
effective-theory side have to be taken into ac- 
count. On the SM side we have to consider of 
the order of 1000 three-loop diagrams. One 
of which is shown in Fig. Obviously, when 
the virtual top quark is present in the open 
fermion line, we have to deal with three-loop 
vacuum integrals involving two mass scales, 
m t and M\y. Also in the charm-quark sector 
such two-scale integrals are encountered when 
closed top-quark loops arise on the virtual 
gluon lines. At present, complete three-loop 
algorithms exist for vacuum integrals involv- 
ing only a single mass scale. We have reduced 



our calculation to such integrals by perform- 
ing expansions around the point m t = M\y 
and for m t 3> Mw- In the latter case, the 
method of asymptotic expansions of Feynman 
integrals has been applied [[13 which has been 
implemented in the C++ program exp [I14|. At 
the physical point where M\y /m t ~ 0.5, both 
expansions work reasonably well. This can be 
seen in Fig. [3] where the one-, two- and three- 
loop results for CV, defined through 






Q(0) 



47T 7 



(4) 



(Q = c,t), are shown as functions of y = 
Mw/m t (fJ,o). (The plots for Cs show an ana- 
log behaviour and can be found in Ref. [|BJ.) 

The variable y changes from to 1, i.e. both 
starting points of our expansions are present 
in the figures. Note the relatively narrow 
ranges of the coefficient values on the verti- 
cal axes. The large mt expansions (up to y s ) 
are depicted by the dot-dashed lines, while the 
expansions around m t = Mw (up to (1 — y 2 ) s ) 
are given by the dashed ones. In the one- and 
two-loop cases, solid curves show the exact 
results. The vertical strips mark the experi- 
mental values for y. 

Comparing the three curves in the one- and 
two-loop cases, one can conclude that a com- 
bination of the two expansions at hand gives a 
good determination of the studied coefficients 
in the whole considered range of y. 

Although we do not know the exact curves 
in the three-loop case, the same pattern seems 
to repeat. In fact, the charm-sector expan- 
sions perfectly overlap in the physical region. 
In the top sector, one can (conservatively) 
conclude that 



yt(3) 



C^ o> (n = m t ) = 12.05 ±0.05, 



c l (Mo = m t ) 



-1.2 ±0.1, 



(5) 
(G) 



which is perfectly accurate for any phe- 
nomenological application. Let us note that a 
change of C 7 (no = rn t ) from 12 to 13 would 
affect the b — > sj decay width by only 0.02%, 
while a similar variation of C 8 (/io = m t ) 
would cause even a smaller effect. More de- 
tailed results can be found in Ref. [|Hj- 

3. Fermionic corrections to Pi, P 2 , Pj 

and P& 

The fermionic NNLO corrections of 
0(a 2 n f ) to Pi, P 2 , P 7 and P 8 (cf. Ref. [ 



|H] for their definition) require the evaluation 
of two- and three-loop diagrams where a light 
fermion loop is inserted into the gluon line. 
Sample diagrams are shown in Fig. 0] 

The computation of the virtual corrections 
proceeds along the same lines as at NLO [ 
I15| . The essential difference is a modified 
gluon propagator which gets an additional e- 
dependent contribution l/p 2e arising from the 
light fermion loop. 

The amplitudes of the contributing dia- 
grams are constructed from one-loop build- 
ing blocks describing the light-fermion inser- 
tion in the gluon propagator and the charm 
quark mass-dependent 6sg7-vertex. After- 
wards Feynman parameters are introduced. 
In the case of P7 and P$ the correspond- 
ing integrations can be performed analytically 
which is not possible for Pi and P 2 . For the 
latter it is convenient to switch for the prop- 
agator involving the charm and the bottom 
quark mass to the Mellin-Barnes representa- 
tions. This leads to a natural expansion in 
m c /m,b which in practice converges very fast; 
the term of 0((m 2 /m 2 ) 3 ) already leads to a 
negligible contribution. 

The practical calculation in the case of Pj 
is slightly more involved as only the combina- 
tion of the virtual corrections with the gluon 
brcmsstrahlung and the quark-pair emission 
process is infrared safe. We decided to reg- 
ulate the infrared singularities in the indi- 
vidual pieces by introducing a finite strange 
quark mass, m s , and a mass of the quark in 
the fermion bubble, rrif. This enables several 
checks on intermediate results. In particular 
we could show that the sum of the virtual and 
the gluon bremsstrahlung corrections are fi- 
nite in four dimensions in the limit m s — > 
and for fixed rrif. The remaining dependence 
on rrif is canceled after including also the con- 
tribution from the quark-pair emission. Ana- 
lytical results for the individual pieces can be 
found in Ref. [ Ell , even for the case where 
a cut-off in the photon energy is introduced 
in the final phase space integration, which is 
important for the comparison with the exper- 
iment. 

The numerical effect of the 0(a 2 nf) cor- 
rections is shown in Fig. OD The correc- 
tions are moderate and amount to —3.9% for 
H = 3.0 GeV and to +3.4% for \x = 9.6 GeV. 
Fig.[S]iriustrates that the /x dependence of the 
0(a 2 nf) improved prediction for the branch- 
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Figure 3. The coefficients C-f (l^o) as functions of y = M\\r/mt(no)- The (blue) dot-dashed lines 
correspond to their expansions in y up to y 8 . The (red) dashed lines describe the expansions in 
(1 — y 2 ) up to (1 — y 2 ) 8 . The (black) solid lines in the one- and two-loop cases correspond to the 
known exact expressions. The (yellow) vertical strips indicate the experimental range for y. 



ing ratio is somewhat flatter than in the NLO 
case if we restrict ourselves to jx > 4 GeV. 
This is a welcome feature of our result, how- 
ever, in general we cannot expect to reduce 
the n dependence as the solid curve only rep- 
resents a part of the NNLO result. Indeed, we 
obtain a stronger /x-dependence in the region 
below 4 GeV. 

4. Summary 

In this contribution we reported on the 
first NNLO calculations to the branching ra- 
tio BR(B — > X s j). In particular, the three- 
loop matching calculations for CV and C& have 
been described [ |Sj ■ Together with the two- 
loop calculation for the other relevant opera- 
tors calculated several years ago [|5] this com- 



pletes the first step towards a complete NNLO 
analysis. 

As far as the corrections to the operator ma- 
trix elements are concerned up to now only 
the fermionic corrections to the numerically 
most important operators Pi (i — 1,2,7,8) 
have been evaluated [EIj ■ The remaining (vir- 
tual and real) corrections are not yet avail- 
able. The same is true for the evaluation of 
the anomalous dimension to NNLO where di- 
agrams up to four loops, from which the pole 
parts are needed, have to be computed. 

It should be stressed that only the complete 
NNLO calculation, outlined after Eq. Q, 
leads to a scheme and scale independent result 
and is thus able to reduce the charm quark 
mass ambiguity. 







Figure 4. Sample diagrams contributing to the matrix elements of the operators Pi (i — 1, 2, 7, 
In the case of Pr next to the virtual corrections also the real radiation has to be considered. 
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Figure 5. The branching ratio BR(_B — » A s 7) as a function of the renormalization scale \i where 
the plot on the right is an enlargement of the one on the left. The dash-dotted curve represents the 
NLO approximation and the solid curve includes the corrections of 0(p? s nf). For illustration in 
the left plot the latter are also shown for the case where the corrections to P\/2 (Pt) are set to zero 
which corresponds to short-dashed (long-dashed) curve. A photon energy cut of E cu t = rrib/20 is 
used, which corresponds to 8 = 0.9. 
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